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Radial-velocity measurements and sine-curve fits to the orbital radial velocity 
variations are presented for ten close binary systems: DU Boo, ET Boo, TX Cnc, 
V1073 Cyg, HL Dra, AK Her, VW LMi, V566 Oph, TV UMi and AG Vir. By 
this contribution, the DDO program has reached the point of 100 published radial 
velocity orbits. The radial velocities have been determined using an improved 
fitting technique which uses rotational profiles to approximate individual peaks 
in broadening functions. 

Three systems, ET Boo, VW LMi and TV UMi, were found to be quadruple 
while AG Vir appears to be a spectroscopic triple. ET Boo, a member of a 
close visual binary with P vis = 113 years, was previously known to be a multiple 
system, but we show that the second component is actually a close, non-eclipsing 
binary. The new observations enabled us to determine the spectroscopic orbits 
of the companion, non-eclipsing pairs in ET Boo and VW LMi. The particularly 
interesting case is VW LMi, where the period of the mutual revolution of the two 
spectroscopic binaries is only 355 days. 

While most of the studied eclipsing pairs are contact binaries, ET Boo is com- 
posed of two double-lined detached binaries and HL Dra is single-lined detached 
or semi-detached system. Five systems of this group were observed spectroscopi- 
cally before: TX Cnc, V1073 Cyg, AK Her (as a single-lined binary), V566 Oph, 
AG Vir, but our new data are of much higher quality than the previous studies. 

Subject headings: stars: close binaries - stars: eclipsing binaries - stars: variable 
stars 



1. INTRODUCTION 

This paper is a continuation of a series of papers (Papers I - VI and VIII - X) of radial- 
velocity studies of close binary stars (Lu & Rucinski 1999; Rucinski & Lu 1999; Rucinski et 
al. 2000; Lu et al. 2001; Rucinski et al. 2001, 2002; Rucinski et al. 2003; Pych et al. 2004; 
Rucinski et al. 2005) and presents data for the tenth group of ten close binary stars observed 
at the David Dunlap Observatory. For technical details and conventions, and for preliminary 
estimates of uncertainties, see the interim summary paper Rucinski (2002a, Paper VII). 

In this paper, we make use of broadening functions extracted not only from the region of 
the Mg I triplet at 5184 A, as in previous papers, but also from two regions containing telluric 
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features centered at 6290 A and 6400 A. These experimental setups were used because of 
concerns about flexure effects in our spectrograph. While this experiment provided a good 
check on the stability of our radial-velocity system and - to a large extent - alleviated our 
concerns, we found that the stellar lines in these two regions were generally too weak to 
replace the 51 84 A feature on a routine basis. The broadening functions (from now on called 
BF's) based on the 6290 A and 6400 A observations were rather poorly defined, especially for 
earlier spectral types; this was mostly due to the low efficiency of our diffraction grating in 
the red region. As the result, the BF's for AG Vir and DU Boo were poor, with the secondary 
component almost undetectable. Thus, in the end, we have returned to the 5184 A region for 
the subsequent observations. The flexure tests based on our telluric-lines template (Regulus 
= HD87901) have shown that the standard wavelength calibrations provide a reasonable 
stability of our radial- velocity system with the largest deviations within ±3 km s" 1 . The 
broadening functions obtained in the red region were used in the present study only to 
augment the data for the quadruple systems, ET Boo, VW LMi and TV UMi, and only for 
observations at critical orbital phases of long-period systems when any spectrum was of use 
in defining the orbit. 

In August 2005, a new grating with 2160 1/mm was acquired to replace the previously 
most frequently used grating with 1800 1/mm which after many years of use lost its original 
efficiency. Thus, unfortunately, due to the changes described above, in combination with 
the necessarily very extended time coverage for triple and quadruple systems, the present 
dataset is the least homogeneous since the start of this series of studies. The BF's used here 
were extracted from spectra obtained with four different CCD detectors and two different 
diffraction gratings. This lack of homogeneity does not seem to affect the final data which 
have uncertainties similar to previously reported in this series of investigations. 

Selection of the targets in our program remains quasi-random: At a given time, we 
observe a few dozen close binary systems with periods usually shorter than one day, brighter 
than 10-11 magnitude and with declinations > —20°; we publish the results in groups of 
ten systems as soon as reasonable orbital elements are obtained from measurements evenly 
distributed in orbital phases. Whenever possible, we estimate spectral types of the program 
stars using our classification spectra obtained with a grating of 600 1/mm over a range of 
635 A or 890 A (depending on the CCD detector) centered at 4200 A. Our classifications 
are based on comparison with several spectral standards of the MK system observed on the 
same nights. They are compared with the mean (B — V) color indices usually taken from 
the Tycho-2 catalog (H0g et al. 2000) and the photometric estimates of the spectral types 
using the relations published by Bessell (1979). 

The radial velocity (hereafter RV) observations reported in this paper have been col- 
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lected between June 1997 and September 2005. The ranges of dates for individual systems 
can be found in Table 1. The present group contains 3 quadruple systems, ET Boo, VW LMi 
and TV UMi, whose complex nature had been noticed several years ago, but whose full or- 
bital solutions required extended monitoring. 

This paper is structured in a way similar to that of previous papers, in that most of the 
data for the observed binaries are in two tables consisting of the RV measurements in Table 1 
and their preliminary sine-curve solutions in Table 2. RVs and corresponding spectroscopic 
orbits for all ten systems are shown in Figures 1 to 3. In this paper we changed the way how 
RVs are determined from the broadening functions: Instead of Gaussians, we now use single 
or double rotational profiles to fit the peaks in the broadening functions. This approach, 
which is described in Section 3, gives much better results with smaller random errors. We 
stress that this is still not a full modeling of the broadening function shape (as attempted in 
Rucinski (1992); Lu & Rucinski (1993); Rucinski et al. (1993)) which would be an optimal 
approach, but a convenient and better working (than Gaussians) tool. The measured RVs 
are listed in Table 1 together with weights, determined as 1/cr 2 , as based on individual 
determinations of central component velocity. This weighting scheme, which accounts for 
differences in relative quality of observations, improves the quality of the orbital solutions. 

The data in Table 2 are organized in the same manner as in previous papers. In addition 
to the parameters of spectroscopic orbits, the table provides information about the relation 
between the spectroscopically observed epoch of the primary-eclipse T and the recent pho- 
tometric determinations in the form of the O — C deviations for the number of elapsed 
periods E. For HL Dra the reference ephemeris is taken from the Hipparcos Catalogue; for 
DU Boo from Pribulla et al. (2005); for the rest of the systems, ephemerides given in the 
on-line version of "An Atlas O-C diagrams of eclipsing binary stars" 4 (Kreiner 2004) were 
adopted. Because the on-line ephemerides are frequently updated, we give those used for 
the computation of the O — C residuals below Table 2 (status as in February 2006). 

The values of T refer to the deeper eclipse which for W-type systems corresponds to the 
lower conjunction of the more massive component; in such cases the epoch is a non- integer 
number. In the cases of ET Boo and VW LMi, where observations covered several years and 
photometric data have been rather scanty, we optimized not only T , but also the orbital 
period. 

Table 2 contains our new spectral classifications of the program objects. Independent 
classification was done for all systems except TX Cnc. Section 2 of the paper contains brief 
summaries of previous studies for individual systems and comments on the new data. The 
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novel technique of fitting the rotational profiles to peaks in the BF's is described in Section 3. 
Examples of BF's of individual systems extracted from spectra observed close to quadratures 
are shown in Fig. 4. 

Similarly as in our previous papers dealing with multiple systems, RVs for the eclipsing 
pair were obtained from BF's with the additional peaks removed. This task was performed 
by first fitting multiple Gaussian profiles to the combined BF's and then removing the 
signatures of the third (and sometimes fourth) component. While the final RVs of the close 
pair were determined by rotational profile fitting to such "cleaned" profiles, the velocities of 
well separated and slowly rotating components of the additional components were determined 
by the Gaussian fits (Table 3). Because the BF technique actually produces Gaussians for 
intrinsically sharp signatures with er ~ 15 km s -1 , this approach is internally consistent. 



2. RESULTS FOR INDIVIDUAL SYSTEMS 
2.1. DU Boo 

Photometric variability of DU Boo was discovered by the Hipparcos satellite (ESA 
1997) where the star was classified as ellipsoidal variable of the A2 spectral type. Later, 
Gomez- Forrellad & Garcia-Melendo (1997) observed the system photometrically and found 
that it is an eclipsing binary with a large O'Connell effect amounting to the difference 
Max. 1 1 — Max. I = 0.10 mag in the V passband. It is interesting to note that light-curve 
asymmetry and the associated surface inhomogeneities have been very stable since the time 
of the Hipparcos discovery; this indicates that solar-type dark-spots paradigm does not 
apply in this case. Recently, Pribulla (2004) analyzed the UBV photometry and found that 
DU Boo is a relatively long-period (1.0559 day) contact binary showing total eclipses; the 
derived photometric mass ratio was found to be q = 0.194(2). Our spectroscopic mass ratio 
q = 0.234(35) is consistent with the photometric determination, which documents reliability 
of photometric mass ratios derived from timing of the inner eclipse contacts for contact 
binaries showing total eclipses (Mochnacki & Doughty 1972a,b). 

The large O'Connell effect is reflected in the extracted broadening functions of DU Boo. 
While the primary component shows undisturbed broadening functions around quadratures, 
close in shape to the theoretical rotational profiles (Section 3), the BF profile for the sec- 
ondary is always very deformed. This causes distortions of the observed RV curve and 
adversely affects the solution for the mass-center velocity. The peaks of the BF's are not 
fully separated supporting the photometric solution of DU Boo as a contact binary. The 
system is clearly of the A-type with the more massive component eclipsed at the deeper 
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minimum. 

By combining our spectroscopic results with the inclination angle % = 81.5° (Pribulla 
2004), we obtain the total mass of the system M 1 + M 2 = 2.56 ± O.O7M . Our new spec- 
tral type estimate of A7V is definitely later than the spectral type given in the Hipparcos 
catalogue (A2). The mean Tycho-2 color (B — V) = 0.31 is in better accord with our deter- 
mination of the spectral type and requires only a small interstellar extinction. The orbital 
period of the system (1.0559 days) is rather long for a contact binary of A7V spectral type 
which indicates that the components of DU Boo may be evolved. The Hipparcos parallax 
is relatively small it = 2.58 ± 1.03 and not precise enough for determination of the system 
luminosity. 

2.2. ET Boo 

The photometric variability of the star was discovered by the Hipparcos satellite (ESA 
1997); where it is cataloged as a (3 Lyrae eclipsing binary of the F8 spectral type. ET Boo is 
a known member of the visual pair COU 1760 (ADS 14593+4649), with the orbital period 
about 113 years and the magnitude difference of AV^ = 0.86 (see Sixth Catalog of Orbits of 
Visual Binary Stars, currently available only in the electronic form 5 ). 

The close, eclipsing binary producing the light variations (from now on: stars 1 and 2) is 
the brighter component of the visual pair. The observed separation of the visual components 
was 0.1 — 0.2 arcsec during the available astrometric observations (Mason et al. 2001) which is 
much less than the typical seeing at the DDO of 1 - 3 arcsec; therefore, the spectra of both 
components were observed simultaneously. We discovered that the broadening functions 
(Figures 4 and 5) show an occasional splitting of the third-component peak, indicating that 
it is in fact a close, non-eclipsing binary (from now on: stars 3 and 4) with a very strongly 
elliptical orbit. Thus, the system is a hierarchical quadruple with both components of the 
visual 113 year period system being double-line (SB2) close binaries. Our estimates of the 
combined brightness of the third and fourth components are L 34 / (Li + L 2 ) = 0.35 ± 0.02 for 
the spectral region around 5184 A and L 34 /(Li + L 2 ) = 0.33 ± 0.02 at 6400 A during the 
maximum light of the closer pair. 

The RV data for the close binary were handled in the standard way, by first removing the 
peaks (by preliminary fitting of Gaussian profiles) of the second binary and then measuring 
the positions of the RV peaks for the close pair. The novelty here is that instead of Gaussians, 
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as in previous papers, we used the double rotational profiles for the close pair (Section 3). 
The results indicate that the brighter component of ET Boo is a semi-detached or more 
likely a detached binary with a relatively large mass ratio, q = 0.884. 

Only a small fraction of the available BF's show splitting of the visual companion 
peaks (components 3 and 4); this property gave us a hint of a highly eccentric orbit, but 
also crucially helped in finding the orbital period. Because these stars have very similar 
brightness, the period ambiguity was resolved by consideration of the RV differences between 
the components. A preliminary orbital period found by trigonometric polynomial fitting to 
the data was later refined by the spectroscopic orbit solution (Table 4, Fig. 6) to give 31.521 
days. The systemic velocity of the second binary Vq = —24.15 ± 0.44 km s -1 is close to 
the systemic velocity of the close pair V = —23.52 ± 0.52 km s" 1 confirming the physical 
association of the two binaries. Since the orbital period of the wide pair is 113 years, no 
orbital motion can be expected to be detected in the 5 years of our observations. 



2.3. TX Cnc 

The photometric variability of TX Cnc, an apparent member of the Praesepe open 
cluster, was first announced by Haffner (1937). Whelan et al. (1973) were the first to 
obtain good simultaneous fit to the photometric and spectroscopic data of the system on 
the assumption of the Roche model. A preliminary analysis of the DDO observations was 
published in a PhD thesis (Blake 2002). The RVs used there have been later re-determined 
using the rotational profile fitting. Surprisingly, the broadening functions show a shape of a 
practically detached binary (see Fig 4), although the orbital velocities are quite typical for 
a contact, W-type system. 

The spectroscopic elements of TX Cnc determined by Whelan et al. (1973) were Vq = 
26.6 ± 3 km s _1 , K 1 = 117.3±3 km s _1 , K 2 = 189.8 ± 4 km s _1 , giving q = 0.62. A later RV 
solution published by McLean & Hilditch (1983), V = 29 ± 6 km s -1 , K 1 = 96 ± 8 km s" 1 
and K 2 = 181 ± 11 km s" 1 , was based on only 8 photographic spectra. Our spectroscopic 
elements suggest a smaller mass ratio q = 0.455 ± 0.011 than derived before and a larger 
total mass, (Mi + M 2 ) sin 3 i = 1.330 ± O.O12M ; both changes are rather characteristic for 
an improved quality of the spectroscopic observations. 

TX Cnc is of particular interest for our understanding of the evolution of contact binaries 
because it appears to belong to Praesepe, which is one of the youngest (900 Myr) open 
clusters containing such binaries (Rucinski 1998). All indications point at an advanced age 
of contact binaries so that confirmation of the membership of TX Cnc to Praesepe may 
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provide a much needed lower limit on the time needed to form such binaries. Unfortunately, 
the system's parallax was not measured during the Hipparcos mission so that the membership 
must be judged by less direct means. Our radial velocity data giving Vq = 34.0 ±0.5 km s -1 , 
are fully consistent with the mean velocity of the Praesepe cluster, V = 34.53 ± 0.12, and 
the velocity dispersion of its spectroscopic binaries, a = 0.40 km s _1 (Mermilliod & Mayor 
1999). 

In the PPM Catalogue, Roser & Bastian (1988) assign TX Cnc a parallax of 5.21 ±0.79 
mas. Hipparcos astrometry of Praesepe analyzed by van Leeuwen (1999) gives a cluster 
parallax of 5.32 ± 0.37 mas, i.e. a distance modulus of m — M = 6.37 ± 0.15. With the 
maximum brightness of TX Cnc of V max = 10.0, we get My = 3.63 ± 0.15 which is in perfect 
agreement with the absolute magnitude estimated from the Rucinski & Duerbeck (1997) 
calibration giving My = 3.60, for assumed (B — V) = 0.54 corresponding to spectral type 
F8V. 

This excellent agreement in radial velocity, parallax and luminosity distance is supported 
by proper motion data. The careful photographic study by Jones & Stauffer (1991) showed 
that TX Cnc has a proper motion of (fj, x = — 4.4±3.2, fi y = 0.5±1.4) mas yr" 1 relative to the 
centre of motion of Praesepe, leading Jones & Stauffer (1991) to assign a 99% probability 
that TX Cnc belongs to Praesepe. The velocity dispersion contributes less than 1 mas yr -1 . 
The Tycho-2 work of H0g et al. (2000) yields (/i a cos5 = -36.2 ± 1.2, = -11.2 ±1.3) mas 
y r -i f Qr aDS0 i u ^ e proper motion of TX Cnc, compared with the mean centre of mass 
motion of (— 35.7±0.4, — 12.7±0.3) mas yr -1 found by van Leeuwen (1999). We can therefore 
be very confident that TX Cnc is a member of Praesepe, and hence it is an important system 
for testing theories of contact binary formation and evolution. 

2.4. V1073 Cyg 

The bright A-type contact binary V1073 Cyg was a subject of several photometric stud- 
ies (Sezer 1996; Morris & Naftilan 2000). There exist also two spectroscopic investigations: 
FitzGerald (1964) published a spectroscopic orbit for the primary component with a marginal 
detection of the secondary, with V = — 8 ± 3 km s _1 , K\ = 66 ± 4 km s _1 and the mass 
ratio of q = 0.34. The author found a small eccentricity of the orbit, e = 0.115 ± 0.053, 
which cannot be significant when one applies arguments of Lucy & Sweeney (1971, 1973) on 
statistics of eccentricity determinations. 

Ahn et al. (1992) obtained an apparently more reliable, circular orbit solution for both 
components, with the resulting spectroscopic elements, V = —0.8 ±1.1 km s -1 , K x = 
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66.7 ± 1.3 km s" 1 and K 2 = 210.2 ± 1.2 km s" 1 . Our results, V = -6.85 ± 0.50 km s _1 , 
Ki = 65.53 ± 0.64 km s -1 and K 2 = 218.9 ±1.5 km s -1 , are superior to the previous 
ones thanks to the BF extraction technique and the rotational profile fitting. The relatively 
large formal rms errors are mainly caused by the simple sine curve solution used by us. 
When combined with a high-precision light curve, the BF modeling can provide high-quality 
absolute parameters of the system. 

Our spectral type estimate, FOV, is much later than original classification of FitzGerald 
(1964) who estimated the spectral type as A3Vm; it confirms the classification of Hill et al. 
(1975). The Tycho-2 color B — V = 0.375 and our spectral type indicate a non-negligible 
amount of reddening. The Hipparcos parallax, 7r = 5.44 ±0.95 mas is not sufficiently precise 
to draw conclusions on the physical properties of the system. 

2.5. HL Dra 

Variability of HL Dra was detected during the Hipparcos mission. The system was 
classified as a (3 Lyrae eclipsing binary with the orbital period of 0.944276 days. The primary 
component is of the A5 spectral type. No ground-based photometric study of the system has 
been published yet. Also no recent minima after the Hipparcos mission are available. Our 
time of the spectroscopic conjunction shows only a small shift (±0.0047 days) with respect 
to the time predicted by the original Hipparcos ephemeris so that the orbital period of the 
system appears to be very stable. 

We have not been able to detect spectral signatures of the secondary component in our 
data. The broadening functions close to quadratures show only small humps on both sides 
of the primary peak which cannot be identified with the secondary component because they 
do not show any orbital motion. The system is clearly a detached or semi-detached pair with 
a low luminosity secondary component. 

HL Dra was observed during two seasons in two different wavelength regions, in 2004 at 
5184 A and in 2005 at 6290 A. The latter dataset is of relatively poor quality due to the low 
number and weakness of spectral lines in the red spectral region of an early type star. The 
orbital single-line solutions resulting from the two datasets are in a good accord except for 
the center of mass velocity of Vq = — 29.3±0.4 km s _1 for the 2004 data and Vq — — 36.5±0.6 
km s _1 for the 2005 data. The shift is well outside the formal errors and may be caused by a 
motion of the eclipsing pair around the barycenter with a third body. The 2004 data are of 
much better quality so the orbital parameters listed in Table 1 correspond to this dataset. 

Our new spectral type determination is slightly later, A6V, than previously published. 
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The Tycho-2 color index (B — V) — 0.222 corresponds to the A8V spectral type so that the 
reddening appears to be small. 

2.6. AK Her 

AK Her is W UMa-type contact binary discovered by Metcalf (see Pickering (1917)). 
It is the brighter component in the visual binary ADS 10408. The companion, located at 
a separation of 4.2 arcsec, is 3.5 mag fainter than AK Her at its maximum light. The 
position angle 323° is almost perpendicular to our fixed, E-W spectrograph slit so that this 
component was not detectable in the broadening functions. 

The system is known to show a cyclic variation in the moments of eclipses which is 
probably caused by the light-time effect induced by an undetected companion on an orbit 
of about 57 years (Awadalla et al. 2004). The perturbing star cannot be identified with the 
known visual companion and must be much closer to the binary. The complex multiplicity of 
the system is supported by the Hipparcos astrometry by the following: (i) The system shows 
a stochastic astrometric solution (the X flag in the catalog field H59); (ii) It is suspected 
not to be single (the S flag in H61), (iii) The trigonometric parallax of 10.47 ± 2.77 mas 
has a much too large error for the brightness of the system. Our individual spectroscopic 
data do not show any contribution from this putative third (or rather fourth) component. 
It is possible that such a companion will be seen in a detailed analysis of averaged spectra 
(D'Angelo et al. 2006), but this approach is outside the scope of the present paper. 

Our RV solution is the first to treat the star as a double-lined binary system (SB2). 
Sanford (1934) observed the RV curve of the primary component and determined f(m) = 
0.0208 M . His RV solution (V = -13 km s _1 , K 1 = 79 km s" 1 ) is fairly consistent with 
our solution. 

Our spectral classification gives an earlier spectral type for the system, F4V, than previ- 
ously discussed, F8V. It is not fully consistent with the Tycho-2 color index (B — V) = 0.490 
and implies some interstellar reddening. 

2.7. VW LMi 

The photometric variability of VW LMi was found by the Hipparcos mission. It was 
classified in the Hipparcos catalogue (ESA 1997) as a W UMa type eclipsing binary with an 
orbital period 0.477547 days. The first photometric observations of the system were published 
by Dumitrescu (2000). Later the light curve of the system was analyzed by Dumitrescu 
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(2003), who found the mass ratio q ph = 0.395 and inclination i = 72.4°; we show later that 
these values are incorrect as they do not take into account the presence of the relatively 
bright binary companion. 

VW LMi has been observed spectroscopically at DDO since 1998. It was realized from 
the beginning that the system is a quadruple one, consisting of two spectroscopic binaries. 
While the eclipsing pair can be identified with the short-period contact binary, the second 
spectroscopic binary is a detached one with the period of about 7.9 days (see below how we 
arrived at the more exact value). The light contribution of the second spectroscopic binary 
at the maximum brightness of the contact pair is L^/{L\ + L<i) = 0.42. 

The study of quadruple system VW LMi is complicated by the mutual orbital motion 
of both binaries so that changes in the respective Vq values cannot be neglected. Another 
complication is the similar brightness of components of the second, non-eclipsing binary 
making derivation of its orbital period very difficult. We worked first with RV differences 
of the components of the second pair to find its orbital period. Trigonometric polynomial 
fits to the data led to the orbital period P 34 = 7.9305 days, which explained all data very 
well. An attempt to find the orbit of the contact pair (after removing the contribution of 
the second binary) resulted in poor quality of the spectroscopic orbit. In fact, the residuals 
from both preliminary orbits showed a clear anti-correlation between the velocities of the 
contact pair (components 1 and 2) and of the non-eclipsing binary (components 3 and 4) 
indicating the mutual orbital motion of the two binaries. The period analysis revealed only 
one feasible period of about 355 days, hence close to one year. Since the data span 7 years 
and the observing season for VW LMi is from late November to mid May, the phase coverage 
of the mutual orbit is partial and has gaps. 

A further improvement of all three orbits was achieved by simultaneous fits to all four 
datasets of the RVs of the form: 

Vi = V + (-lfKilej cosujj + cos(^ + Vj)] + (-l) j 'K 2j -i, 2 j[e3 coscj 3 + cos(cj 3 + i/ 3 )] 

where Vq is the mass-center velocity of whole quadruple system, are the respective velocity 
semi-amplitudes of individual components, is orbital eccentricity, ujj is longitude of the 
periastron and Uj is the true anomaly. The index % corresponds to the component number 
{% — 1 — 4), while the index j takes the value of 1 for the contact binary, 2 for the detached 
binary and 3 for the mutual orbit of the two systems. Thus, for the components of the 
contact binary, j = 1 and i — 1,2 while for the components of the detached binary, j = 2 
and % = 3,4. T^j-i^j for j—1 should be read as K\ 2 while for j = 2 as K34 where K12 and 
K34 are the semi-amplitudes of the RV changes of the mass centers of the contact and the 
detached binary, respectively. 
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All results of the simultaneous fits are presented in Table 5 while the orbital elements 
of the contact pair are also given with the remaining binaries of this study in Table 2. 
For simplicity, all measurements were assigned the same weight although the velocities for 
the detached pair were determined by the Gaussian fits while those of contact pair by the 
rotational profile fits. The sine curve fits to the data for the contact binary, corrected for the 
motion on the outer orbit, are shown in Fig. 2. While the secondary component is usually 
not blended with the peaks of the second spectroscopic binary, the primary of the contact 
pair is always visible projected against the "background" of the profiles of the 3rd and 4th 
component. This circumstance caused an enhanced scatter of the velocities of the primary 
component. 

The corrected RVs of the second binary with the corresponding fits are plotted in Fig. 7. 
The final orbital period for this binary is 7.93044 days and the orbit is nearly circular. The 
velocities of all four components corrected for the corresponding orbital motions in the inner 
orbits and their best fits are plotted in Fig. 8. These residual RVs represent the orbital 
motion of the mass centers of both binaries. 

Because the outer orbit has a relatively short period of 355 days, it is of interest to inquire 
into the mutual orientation of the three orbits. This can be estimated from projected masses 
of the components in a sort of a "bootstrap" process started with the derived inclination 
of the contact, eclipsing system. A preliminary solution of unpublished photoelectric data 
obtained at the Stara Lesna Observatory of the Astronomical Institute of the Slovak Academy 
of Sciences were used to estimate inclination angle of the eclipsing pair. Fixing the third 
light at L 34 /(Li + L 2 ) = 0.42 (see above) and using the spectroscopic mass ratio of q = 0.42 
led to the inclination angle %vi = 80.1° ±0.2°. This is, as expected, a much larger inclination 
than the one obtained by Dumitrescu (2003) (i = 72.4°) without an assumption of the third 
light. Using our estimate of the inclination angle and the projected total mass of the contact 
pair (Mi + M 2 ) sinii 2 = 2.28 M , we obtain {M l + M 2 ) = 2.39 M . The outer, 355 day orbit 
defines the mass ratio for the two pairs, (Mi + M 2 )/(M 3 + M4) = 1.09. Therefore, the true 
(not the projected) mass of the second spectroscopic binary is (M3 + M4) = 2.19 M . Using 
the projected mass (M 3 + M 4 ) sini 34 = 1.79 ± 0.03 M , we estimate the inclination of the 
orbit of the second pair to be about 69°. The outer, 355 day orbit is even less inclined to 
the sky since with (Mi + M 2 + M 3 + M 4 ) = 4.58 M and the projected total mass of only 
2.67 M , so that we obtain in-u = 57°. Obviously, we could not determine if these values 
are all in the same sense or are complements to 180°; a determination of the sense of the 
revolution could only come from interfermetric observations. 

The non-eclipsing, detached binary, with (M 3 + M 4 ) = 2.19 M , is composed of two 
almost identical (g 34 = M 3 /M 4 = K 4 /K 3 = 0.980 ± 0.017), most probably main sequence 
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stars. Their masses correspond to about F9V - GOV and their similarity is also reflected in 
the luminosity ratio / 3 // 4 pa 1.04. The evolutionary status of components can be guessed from 
a comparison of their rotational and orbital velocities. If we assume a synchronous rotation 
and take rotational velocities of the components estimated from Gaussian profile fits to be 
about 12 km s _1 and semi-amplitudes of RVs about 60 km s -1 , we see that fractional radii 
of components are r/a < 0.2. With the semi-major axis of the absolute orbit of about 10 
R , their radii are < 2R Q . The similar spectral type of all components of the quadruple 
system of VW LMi has resulted in a very good quality of extracted BF's, as can be judged 
in Figure 4. 

It is interesting to note that multiplicity of VW LMi was not identified astrometrically 
during the Hipparcos mission in spite of the relative proximity of the system at n = 8.04±0.90 
mas. This is probably caused by the orbital period of the two pairs around each other being 
close to one year thus mimicking the parallactic motion. 

Because of the small size of the mutual (355 day) orbit of only 0.62 AU, chances of 
resolving the two astrometric components are rather small even with advanced techniques 
because the expected maximum angular separation will be only 10 mas. The situation is little 
bit more optimistic with the expected light-time effect of the eclipse timing of the contact 
binary as the expected full amplitude is about 0.0074 days which should be relatively easy 
to detect with the observed photometric amplitude of 0.42 mag. 

2.8. V566 Oph 

The W UMa-type binary V566 Oph was discovered by Hoffmeister (1935). The system 
is bright (V max = 7.46) and therefore it was a subject of numreous previous photometric (for 
references, see Twigg (1979)) and spectroscopic observations. An interval of constant light 
observed during the secondary eclipse indicates the A-type. Mochnacki & Doughty (1972b) 
published the first light curve analysis of V566 Oph based on the Roche model. The total 
eclipses permitted to determine reliable geometric elements: the fill-out F = 1.25 ± 0.05 (/ 
= 0.25), mass ratio q = 0.238 ± 0.005 and inclination of i = 80° ± 2°. 

There exist three previous spectroscopic studies of the system. Two of these are based 
on photographic observations (Heard 1965; McLean 1983) and a more recent one is based 
on the Reticon data (Hill et al. 1989). The latter study used direct fitting of the synthetic 
profiles to the CCF functions. The spectroscopic elements obtained in his study, Vo = 
38.5±1.1 km s _1 , K 1 = 72.6±1.5 km s" 1 and K 2 = 272.9±1.3 km s" 1 , are practically within 
the errors of our results, V = 37.3±0.5 km s" 1 , K x = 71.1±0.7 km s" 1 and K 2 = 270.1±1.1 



-14- 



km s _1 . The current improvement of the orbit is mainly due to the use of the BF extraction 
technique and of the rotational-profile fitting. The orbit can be still improved by taking the 
proximity effects into account. Our new determination of the mass ratio, q = 0.263 ± 0.012 
is in a moderately good agreement with photometric mass ratio of Mochnacki & Doughty 
(1972b), confirming the utility of the photometric approach to systems with total eclipses. 

The orbital period of the system is rather unstable. In spite of the possible light-time 
effect orbit found by Pribulla & Rucinski (2006), we did not find any traces of the third 
component in the extracted BFs. 

V566 Oph is a relatively nearby system with a good Hipparcos parallax, n = 13.98±1.11 
mas. The absolute magnitude determined using the calibration of Rucinski & Duerbeck 
(1997), My = 3.07, using [B — V) — 0.406 from the Tycho-2 catalogue (H0g et al. 2000) 
is in good agreement with the visual absolute magnitude determined from the Hipparcos 
parallax, My = 3.19 ± 0.17. Our new spectral classification, F4V, indicates a slightly later 
spectral type than that F2V found by Hill et al. (1989). 

2.9. TV UMi 

TV UMi is another discovery of the Hipparcos mission. The system was classified as a 
/3 Lyrae eclipsing binary with orbital period of 0.415546 days, although the classification was 
obviously complicated by the low amplitude of the light variations of only about 0.08 mag. 
The eclipses are very wide and of almost the same depth. For that reason we suspected that 
the variability is caused by a contact binary of the W UMa-type. 

Prolonged spectroscopic observations of TV UMi at the DDO showed that the system 
is a quadruple one and consists of two spectroscopic binaries. The second close binary in 
TV UMi is almost as bright as the contact pair, but its components are difficult to analyze 
because of the strong eccentricity of the orbit and the very short duration of periastron 
passages when the spectral signatures could be potentially resolved. In fact, the components 
of the second pair could not be separated in most of our spectra; such a separation took place 
on only three occasions. The largest observed separation of the components of the second 
binary on June 8, 2001 of more than 100 km s _1 had to be close to periastron passage. One 
of the BF's from that night is included in Fig. 3. During this event, the stronger component 
had a more negative RV. A period analysis of line separations indicates a 31.2-days orbital 
periodicity of the second pair. The reliable determination of the orbital parameters and 
confirmation of the preliminary orbital period would require more observations, preferably 
with a larger spectroscopic resolution to separate the components even outside the periastron 
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passages. The TV UMi system resembles ET Boo in that its companion binary also has an 
eccentric orbit with a similar orbital period. Our RV observations cover a shorter time base 
and are less numerous than in the case of VW LMi so that we have not been able to find the 
mutual orbital motion of the binaries. During intervals when we observed the narrow blend 
of the peaks of the third and fourth components, the combined RV was about —15 km s -1 , 
slightly less than systemic velocity of the contact binary (—9.7 km s _1 ). This indicates a 
possible slow orbital motion of both pairs. 

The light contribution of the second binary is large, L^j{L\ + L 2 ) = 0.90. As observed 
during the periastron passage on June 8, 2001, its components have slightly unequal bright- 
ness: L 3 /(Li + L 2 ) = 0.58 and L 4: /(L 1 + L 2 ) = 0.32. Using the measured RVs of the third 
and fourth components during the periastron passage, RV 3 = —67.94 ± 1.13 km s _1 and 
RV4 = 37.21 ± 1.03 km s _1 , and the mean RV of the blend of the two components giving the 
approximate systemic velocity of the second pair, —15 km s -1 , we see that the mass ratio 
of the second binary is close to unity. The observed photometric amplitude of the contact 
system, Am = 0.083, when corrected for light contribution of the companion binary, remains 
small at Am = 0.163. 

It is interesting to note that the system was not detected nor even suspected as a 
multiple one from the Hipparcos astrometric data. Chances to resolve both binaries by direct 
imaging are higher than in VW LMi because no rapid orbital motion of the contact binary 
was observed indicating a longer orbital period and thus a larger separation of both binaries. 
According to the Hipparcos astrometry TV UMi is relatively nearby, with 71 = 7.64 ± 0.78 
mas. 

Clearly, the current data for the second close binary system in TV UMi are inadequate 
for determination of full orbital parameters for the whole system. Such a determination 
would require a long-term monitoring program with one or two spectra obtained per night 
over a period of few months. 

2.10. AG Vir 

Variability of AG Vir was discovered by Guthnick & Prager (1929). Since then it was 
a subject to several photometric investigations (for references, see Bell et al. (1990)). The 
system is very similar to DU Boo (Section 2.1) in that the first light maximum is always 
brighter of the two by about 0.08 mag. A simultaneous photometric and spectroscopic 
analysis of Bell et al. (1990) showed that a reliable determination of the geometrical elements 
was complicated by the strong asymmetry of the light curve. The system is of the A-type 
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with an undeterminable degree of contact (could be marginal or deep). The spectral type is 
A7V (our new determination suggests A5V), so that the observed photospheric brightness 
inhomogeneities may be quite different from the solar-type dark spots. 

The period analysis of Blanco & Catalano (1970) indicated the presence of a light-time 
effect in the eclipse timing caused by a third component on a 40 years orbit. However, later 
observations did not confirm any cyclic behavior, although the observed times of minima 
do show a large scatter. This can be interpreted either by the light-curve variations caused 
by the presence of surface spots or by a light-time effect caused by a body on a short- 
period orbit. In fact, our broadening functions do show a well defined signature of the third 
component with L 3 / (Li±L 2 ) ~ 0.05 (see Fig. 9). However, this star has a rather different RV 
(+2.9 ±2.5 km s~ x ) from the systemic velocity of the contact binary (—10.99 ±0.82 km s _1 ). 
Unfortunately, the uncertainties of the RVs are quite large and only our better-quality 2005 
spectra (as used in this paper) show the presence of this component. A close inspection of 
the published cross-correlation functions of Bell et al. (1990) reveals a marginally defined 
presence of a third component close to the systemic velocity of the binary, but - as expected 
- the definition of the CCF's was inferior to that of the BF's. 

The Hipparcos parallax (1.33 ± 1.18 mas) is too small and inconsistent with the esti- 
mated absolute magnitude for a A5V star, My = ±2.1, and the observed maximum visual 
magnitude, V max = 8.50. It is possible that this inconsistency is caused by a transverse 
motion of the eclipsing pair around the third component. 

3. Measurements of radial velocities using rotational profiles 

While the previous papers of this series (Papers I-VI and VIII-X) reported RV mea- 
surements obtained by Gaussian fits to the extracted broadening functions, in this paper 
we use a novel technique of extracting RVs by rotational profile fits. While neither Gaus- 
sians nor rotational profiles can replace the full modeling of broadening functions which will 
hopefully take place one day, the rotational profiles are equally simple to implement as the 
Gaussians, but offer an improvement in the quality of the RV measurements with a much 
better convergence to the final result and noticeably smaller random errors. 

The broadening function of a rigidly rotating star (Gray 1976) is described by four 
parameters: the overall strength or the amplitude of the BF, a , the central velocity, a± 
(identified with the light centroid velocity of the star), and the half-width ai = Usini; 
an additional parameter is a vertical background displacement, a^, which can be usually 
traced to a different continuum and pseudo-continuum levels for the standard and program 
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stars during the spectral normalization step. For a double-peaked profile, the three first 
parameters (a to a 2 ) appear twice so that the fit involves 7 unknowns. By using auxiliary 
quantities, c± — 1 — u and c 2 = iru/4, the profile can be written as: 



I = a + a>i 



CiJl-(x- ai) 2 /«2 + c 2 (l - (x - ai) 2 /a 2 ) 



/(ci + c 2 ). 



(1) 



The profile depends only slightly on the limb darkening coefficient, u. We assumed 
u = 0.7 in our measurements. As opposed to the Gaussian profile, the rotational profile 
(see an example in Fig. 10) has rather steep sides and - by definition - is exactly zero for 
velocities larger than projected equatorial velocity V rot sini. These properties are crucial 
to the improvement in the determination of the parameter ai, the light centroid for each 
component. 

Although the rotational profiles represent the BF profiles of double-line binaries much 
better than Gaussians, a practical application may encounter the same complications: 



• If the secondary shows a very faint peak in the BF, its width must be usually fixed at a 
reasonable value to improve the stability of the solution and to give consistent results; 

• The BF profiles for contact binaries are rather different from those of single, rigidly 
rotating stars, particularly in the "neck" region between the stars where an additional 
light is present. Obviously, these asymmetric deviations cannot be fitted by rotational 
profiles nor by Gaussians. Normally, this leads to under-estimates of the velocity semi- 
amplitudes. 



In the case of close binary stars, the combined rotational profiles can be applied to broadening 
functions strictly only outside the eclipses as a proper representation of the data should 
involve inclusion of the eclipse and proximity effects (Rucinski 1992). However, a simple 
upper envelope of the two individual profiles works well even during partial eclipse phases. 
We found that the double rotational profiles converge faster to the final results and describe 
the data much better than the Gaussians. This is well illustrated in the case of TX Cnc, 
where a preliminary orbit defined by RVs obtained through a Gaussian fits had almost twice 
as large standard errors of the spectroscopic elements. While rms errors for the velocities 
derived from the Gaussian fits were o\ = 8.4 km s _1 and a 2 = 8.3 km s -1 , the errors for the 
RVs derived by rotational profile fitting are <j\ =4.1 km s" 1 and a 2 — 5.7 km s -1 . 



-18- 



4. SUMMARY 

With the new ten short-period binaries, this paper brings the number of the systems 
studied at the David Dunlap Observatory to a round number of one hundred. The systems 
presented in this paper include three quadruple systems for which we have been collecting 
data for several years in a hope of being able to study variability of RVs on times scales 
ranging from a fraction of the day to several years. This has been achieved for ET Boo and 
VW LMi where we can say much about all components of these hierarchical binaries. ET Boo 
is a known visual binary with the period of 113 years with each component being a close 
binary. VW LMi is a particularly interesting system with the period of mutual revolution 
of both binaries of only 355 days. Starting with a preliminary photometric solution of the 
light curve of VW LMi which gave the orbital inclination of the close binary, we were able to 
determine the orbital inclinations of all involved binaries in this system and thus to derive 
masses of all components. We have been less successful for the quadruple system TV UMi 
where the second pair requires a very prolonged monitoring for analysis of the 31.2 day orbit 
of the second pair. 

We have found that AG Vir appears to be a triple system, although there is inconsis- 
tency in the velocity of its companion. For AK Her, we were able to obtain data free of 
contamination from the known third component; there are indications that this binary has 
another faint companion causing the light-time effect in eclipse timing. The systems DU Boo, 
TX Cnc, V1073 Cyg, V566 Oph are relatively mundane double-lined contact binaries while 
HL Dra is a single-lined binary of an unknown variability type. 

All RVs for close binaries analyzed in this paper have been determined by using a novel 
technique of rotational profile fitting to the broadening functions. This technique, while still 
not perfect in reproducing asymmetries and intricacies of the real BF's, is much advantageous 
and accurate than the Gaussian fitting previously used in our studies. 
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Captions to figures: 

Fig. 1.— Radial velocities of the systems DU Boo, ET Boo, TX Cnc and V1073 Cyg 
are plotted in individual panels versus the orbital phases. The lines give the respective 
circular-orbit (sine-curve) fits to the RVs. ET Boo is quadruple system composed of detached 
eclipsing pair (wide profiles) and detached non-eclipsing pair with a much longer orbital 
period P 31.5 days (narrow peaks). DU Boo, TX Cnc and V1073 Cyg are contact 
binaries. The circles and triangles in this and the next two figures correspond to components 
with velocities V\ and V 2 , as listed in Table 1, respectively. The component eclipsed at the 
minimum corresponding to T (as given in Table 2) is the one which shows negative velocities 
for the phase interval 0.0 — 0.5. Short marks in the lower parts of the panels show phases of 
available observations which were not used in the solutions because of the blending of lines. 
All panels have the same horizontal range, —500 to +500 km s _1 . 

Fig. 2. — The same as for Figure 1, but for HL Dra, AK Her, VW LMi, and V566 Oph. 
HL Dra is a detached single-lined binary while AK Her, and V566 Oph are contact binaries. 
VW LMi is a quadruple systems containing the eclipsing contact binary and a detached 
non-eclipsing binary with a 7.93 days orbit. The study and interpretation of broadening 
functions is complicated by the fact that both binaries revolve in a relatively tight orbit with 
P ks 355 days. We show additional data for VW LMi in Figures 7-8. 

Fig. 3. — The same as for Figures 1 and 2, but for the two remaining systems TV UMi 
and AG Vir. While AG Vir is contact binary containing a third component, TV UMi is 
quadruple system containing a contact binary and a detached non-eclipsing binary with a 
period of about 31.2 days. 

Fig. 4. — The broadening functions (BF's) for all ten systems of this group, selected for 
phases close to 0.25 or 0.75. The phases are marked by numbers in individual panels. For 
the three quadruple systems, ET Boo, VW LMi and TV UMi, we selected the BF's showing 
the second, non-eclipsing pair. While the companion binary of VW LMi has an almost 
circular orbit and its lines are separated during most of the orbit, the lines of the companion 
binaries of ET Boo and TV UMi can be separated only during short intervals during the 
respective periastron passages. 

Fig. 5. — A grayscale plot of broadening functions of ET Boo extracted from spectra taken 
in February, 2005 around 6400 A. The spectra were sorted in time. One bin in X-axis 
corresponds to one spectrum, so the scale is not continuous. The signatures of the close, 
short-period binary are seen as relatively faint features in the background. The second 
binary is visible as a sharp feature in the BF's which split into two components during the 
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periastron passage on the second night of observations, on February 10, 2005. 

Fig. 6. — The radial velocities of the third and fourth component of ET Boo and the 
corresponding fits to an orbital motion with the period of 31.52 days. 

Fig. 7. — The radial velocities of the third and fourth component of VW LMi, plotted in 
a phase diagram with the period of 7.93 days, after being corrected for the motion on the 
outer, 355 day period orbit. 

Fig. 8. — Radial velocities of all four components of VW LMi corrected for orbital motion 
in both close orbits showing motion of the mass centers in the wide, 355 days orbit. The 
phases are counted from the periastron passage. Because the orbital period is close to one 
year, the seasonal observing interval slowly moved in the orbital phases between 1998 and 
2005. 

Fig. 9. — A grayscale plot of the broadening functions of AG Vir. The spectra have been 
sorted in phase, but not rebinned to the proper phase durations. While individual broad- 
ening functions do not show the third component very clearly, the grayscale plot shows a 
continuous, bright ridge line close to the systemic velocity. 

Fig. 10. — A comparison of double Gaussian and rotational fits to the broadening function of 
V566 Oph extracted from spectrum taken close to the orbital quadrature. The RV is better 
defined in the steep portions of the rotational profile branches. 
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Table 1. DDO radial velocity observations of close binary systems (the full table is 

available only in electronic form 



HJD-2,400,000 


Vi 


Wi 


v 2 


W 2 


Phase 




[km s- 1 ] 




[km s -1 ] 






53480.5576 


22.48 


0.35 


-168.48 


0.03 


0.6224 


53480.5803 


25.78 


0.47 


-203.77 


0.03 


0.6439 


53480.6033 


31.53 


1.49 


-209.61 


0.14 


0.6657 


53480.6261 


34.74 


1.09 


-216.95 


0.05 


0.6873 


53480.6490 


33.67 


1.05 


-227.17 


0.08 


0.7090 


53480.6717 


40.69 


1.00 


-233.93 


0.30 


0.7305 


53480.6944 


37.69 


1.19 


-235.97 


0.13 


0.7520 


53480.7172 


38.72 


0.89 


-251.40 


0.13 


0.7736 


53480.7399 


38.37 


1.03 


-237.67 


0.10 


0.7951 


53480.7626 


40.39 


1.32 


-221.31 


0.12 


0.8165 



Note. — The table gives the RVs Vi and associated weights Wi for 
observations described in the paper. The first 10 rows of the table 
for the first program star, DU Boo, are shown. Observations lead- 
ing to entirely inseparable broadening function peaks are given zero 
weight; these observations may be eventually used in more extensive 
modelling of broadening functions. The observations in the clumn V\ 
correspond to the component which was stronger and easier to mea- 
sure in the analysis of the broadening functions; it was not always 
the component eclipsed during the primary minimum at the epoch 
T (see Table 2). The figures should help in identifying which star is 
which. 
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Table 2. Spectroscopic orbital elements 





Tvnp 
rype 


(l+rlPV X~\ Cl TYl OQ 


Vn 


K'-i 




T„ _ o 400 000 


r ^uayby 


Q 




Sp. type 






K 2 


€2 


(O - C)(d) [El 

V v v / V J L J 


(Mi + Mi) sin 3 i 




DU Boo 


EW(A) 


HD 126080 


-13.09(0.86) 


53.59(1.02) 


5.95 


53,494.6896(39) 


1.05588870 


0.234(35) 




A7V 


HIP 70240 




229.32(3.02) 


11.22 


+0.0068 [1173] 


2.477(73) 




ET Boo 


EB 


BD+47 2190 


-23.35(0.51) 


145.61(0.75) 


9.09 


52701.5928(8) 


0.6450398(7) 


0.884(15) 




F7V 


HIP 73346 




164.64(0.93) 


10.09 


-0.0008 [312] 


1.996(17) 




TV Pnr 
_L J\. L'llC 


Hi VV ^ VV ) 




Too.y i i^u.oz ) 


1UU. ( I I^U.oZ ) 


A HQ 


oi, ou i .yoiu^4 j 


U.OO/OOZ 1 o 


U.40D^ll ) 




(F8V) 


BD+19 2068 




221.67(0.91) 


5.66 


+0.0008[-1807.5] 


1.330(12) 




V1073 Cyg 


EW(A) 


HD 204038 


-6.85(0.50) 


66.45(0.61) 


4.81 


53,194.3876(12) 


0.7858506 


0.303(17) 




A9V 


HIP 105739 




219.09(1.40) 


6.63 


+0.0025 [883] 


1.896(25) 




HL Dra 


EB(SBl) 


HD 172022 


-29.36(0.43) 


81.04(0.60) 


1.88 


53,166.8057(12) 


0.944276 






A6V 


HIP 91052 








+0.0047 [4942] 






AK Her 


EW(A) 


HD 155937 


+4.28(0.89) 


70.52(1.12) 


6.90 


53176.3946(19) 


0.4215231 


0.277(24) 




F4V 


HIP 84293 




254.40(2.27) 


13.90 


+0.0009 [1,604] 


1.598(29) 




VW LMi a 


EW(A) 


HD 95660 




105.41(0.83) 


12.68 


51,973.4117(4) 


0.4775499(2) 


0.416(4) 




F5V 


HIP 54003 




253.21(0.84) 


13.58 


-0.0017[-1103] 


2.282(18) 




V566 Oph 


EW(A) 


HD 163611 


-37.33(0.52) 


71.08(0.69) 


3.12 


53,568.6298(4) 


0.4096538 


0.263(12) 




F4V 


HIP 87860 




270.12(1.14) 


6.63 


+0.0037 [2608] 


1.686(17) 




TV UMi 


EW(W) 


HD 133767 


-9.70(0.67) 


116.17(1.04) 


7.57 


52,454.0189(6) 


0.41554935 


0.739(21) 




F8V 


HIP 73474 




157.09(1.19) 


9.98 


+0.0007 [-111.5] 


0.879(12) 




AG Vir 


EW(A) 


HD 104350 


-10.99(0.82) 


93.39(1.06) 


5.38 


53,501.5388(13) 


0.6426507 


0.382(21) 




A5V 


HIP 58605 




244.24(1.97) 


10.81 


+0.0030 [1558] 


2.563(41) 





a VW LMi: The contact binary revolves in relatively short-period orbit around the second pair, therefore systemic velocity Vb is n °t 
defined. See Table 5 for a full description of the system. 



Note. — The spectral types given in column two relate to the combined spectral type of all components in a system; they are given 
in parentheses if taken from the literature, otherwise arc new, The convention of naming the binary components in the table is that the 
more massive star is marked by the subscript "1", so that the mass ratio is defined to be always q < 1. The figures should help identify 
which component is eclipsed at the primary minimum. The standard errors of the circular solutions in the table arc expressed in units 
of last decimal places quoted; they are given in parentheses after each value. The center-of-mass velocities (Vb), the velocity amplitudes 
(Ki) and the standard unit-weight errors of the solutions (e) are all expressed in km s — 1 . The spectroscopically determined moments 
of primary minima are given by To ; the corresponding (O — C) deviations (in days) have been calculated from the available prediction 
on To, as given in the text, using the assumed periods and the number of epochs given by [E], The values of (Mi + M2) sin 3 i are in the 
solar mass units. 

Ephemerides (HJD min - 2,400,000 + period in days) used for the computation of (O-C) residuals: DU Boo: 52256.1254 + 1.0558887; 
ET Boo: 52500.3407 + 0.6450413; TX Cnc: 52500.0407 + 0.38288273; V1073 Cyg: 52500.479 + 0.7858506; HL Dra: 48500.189 
+ 0.944276; AK Her: 52500.2709 + 0.42152292; VW LMi: 52500.1516 + 0.4775505; V566 Oph: 52500.249 + 0.4096538; TV UMi: 
52500.352 + 0.41555; AG Vir: 52500.286 + 0.6426507; 
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Table 3. Radial velocity observations of the third and fourth components of quadruple 
systems (the full table is available only in electronic form) 



HJD-2,400,000 V 3 V 4 

[km s -1 ] [km s -1 ] 



ET Boo 



53455.89155 


-14 


.30 


-42. 


.26 


b 


53455.90414 


-14 


.30 


-43. 


.31 


b 


53455.91510 


-16 


.13 


-42. 


.56 


b 


53455.92738 


-12 


.09 


-36 


,59 


b 


53529.61628 


-89 


.18 


68. 


.01 




53530.73133 


-63 


.86 


44. 


.00 




/W LMi 












50852.78735 


-75 


.68 


51 


,31 




50852.79458 


-75 


.45 


52. 


.17 




50853.69796 


-54 


.16 


29 


,50 




50853.70383 


-55 


.28 


28. 


,85 





Note. — The table gives the RVs Vt for the 
third and fourth components. The typical 10 
rows of the table for the quadruple systems, 
ET Boo and VW LMi, are shown. Observa- 
tions of quadruple systems ET Boo, VW LMi 
and TV UMi leading to entirely inseparable 
broadening function peaks of components of 
the second binary are omitted from the ta- 
ble. Radial velocities determined from par- 
tially blended profiles arc marked by "b" in 
the last column. 
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Table 4. Spectroscopic orbital elements of the second non-eclipsing binary in the 

quadruple system ET Boo 



Parameter 




error 


P34 [days] 


31.5212 


0.0005 


e34 


0.740 


0.011 


uu [rad] 


2.94 


0.03 


T [HJD] 


2,451,354.66 


0.04 


V [km s- 1 ] 


-24.15 


0.44 


K 3 [km s- 1 ] 


40.36 


0.67 


Ki [km s- 1 ] 


57.67 


0.69 


q = K 4 /K 3 


0.70 


0.09 


(as + a 4 ) sini [AU] 


0.191 


0.003 


(M 3 + Mi) sin 3 i [M ] 


0.93 


0.05 



Note. - The table gives spectroscopic ele- 
ments of the second binary in ET Boo: orbital 
period (-P34), eccentricity (634), longitude of the 
periastron passage (w), time of the periastron pas- 
sage (To), systemic velocity (Vo), semi-amplitudes 
of the RV changes (Ks,Ki). Corresponding 
mass ratio q, and projected relative semi-major 
((a 3 + a 4 ) sin i) and total mass ((M 3 + M4) sin 3 i) 
is also given. 
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Table 5. Spectroscopic orbital elements of all three orbits defined by RVs of four resolved 

components of VW LMi 



Parameter 




error 


Contact (eclipsing) pair - circular orbit 




Pro [daysl 


0.47754988 


0.00000020 


T [HJD] 


2,451,973.4117 


0.0004 


K 1 [km s- 1 ] 


105.41 


0.83 


K 2 [km s- 1 ] 


253.21 


0.84 


(Mi +M 2 )sin 3 ii2 [Mq] 


2.282 


0.018 


Detached non-eclipsing pair 






P34 [days] 


7.93062 


0.00014 




0.033 


0.010 


lu 2 [rad] 


5.01 


0.09 


T [HJD] 


2,452,282.44 


0.11 


K 3 [km s- 1 ] 


65.44 


0.76 


K 4 [km s- 1 ] 


64.15 


0.76 


(M 3 + Mi) sin 3 i 34 [M Q ] 


1.788 


0.033 


Mutual wide orbit 






Pi 2 -34 [days] 


355.0 


0.5 


e3 


0.14 


0.03 


w 3 [rad] 


2.35 


0.22 


T [HJD] 


2,452,703 


11 


K 12 [km s- 1 ] 


19.96 


0.84 


^34 [km s- 1 ] 


21.76 


0.79 


(Mi + M 2 + M 3 + M 4 ) sin 3 ii 2 _ 34 [M ] 


2.67 


0.16 


V [km s- 1 ] 


1.29 


0.39 



Note. — The table gives spectroscopic elements for the three observed 
orbits of VW LMi. The designation of parameters is as in the previous 
table. The index "12" refers to the orbit of the contact pair, while the 
index "34" refers to the orbit of the second, detached binary, parameters 
of the mutual orbit of these binaries are indexed as "12-34". The orbit of 
the contact pair is assumed to be circular (ei=0, u>i = n/2). 



